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Abstract. A class of functions involving the divided differences of the psi 
function and the polygamma functions and originating from Kershaw's double 
inequality are proved to be completely monotonic. As applications of these 
results, the monotonicity and convexity of a function involving ratio of two 
gamma functions and originating from establishment of the best upper and 
lower bounds in Kershaw's double inequality are derived, two sharp double 
inequalities involving ratios of double factorials are recovered, the probability 
integral or error function is estimated, a double inequality for ratio of the 
volumes of the unit balls in R"~^ and R" respectively is deduced, and a 
symmetrical upper and lower bounds for the gamma function in terms of the 
psi function is generalized. 



1. Introduction 

Recall [29, Chapter XIII] and [57, Chapter IV] that a function / is said to be 
completely monotonic on an interval / if / has derivatives of all orders on / and 

> (1) 

for X € I and n > 0. The famous Bernstein's Theorem in [57, p. 160, Theorem 12a] 
states that a function / is completely monotonic on [0, oo) if and only if 

/>oo 

fix)= / e~^M/i(s), (2) 
Jo 

where /i is a nonnegative measure on [0, cxd) such that the integral (2) converges for 
all a:; > 0. This expresses that a completely monotonic function / on [0,oo) is a 
Laplace transform of the measure fi. 

Recall also [5, 46, 49] that a positive function / is called logarithmically com- 
pletely monotonic on an interval / if / has derivatives of all orders on / and its 
logarithm In / satisfies 

(-l)^[ln/(x)](^) >0 (3) 

for all G N on /. It was proved explicitly in [10, 42, 45, 46, 53] by different 
approaches that any logarithmically completely monotonic function must be com- 
pletely monotonic, but not conversely. It was pointed out in [10, Theorem 1.1] 
and [24, 50] that the logarithmically completely monotonic functions on [0,oo) are 
those completely monotonic functions on [0, oo) for which the representing measure 
/i in (2) is infinitely divisible in the convolution sense: For each n G N there exists 
a positive measure v on [0, cxd) with n-th convolution power equal to /i. 
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It is well-known that the classical Eulcr gamma function 



T{x) = r 

Jo 



x-\. 



At 



(4) 



for a; > 0, its logarithmic derivative, denoted by i){x) = ^^'^ '^'^ polygamma 

functions V'^'^ (a;) for i e N are several of the most important special functions and 
have much extensive applications in many branches such as statistics, probabil- 
ity, number theory, theory of 0-1 matrices, graph theory, combinatorics, physics, 
engineering, and other mathematical sciences. 

The ratio p|^^^j for .x+p > and x^q > of two gamma fimctions, called Wallis 
function or ratio in the literature, has been investigated since 1948 in [56] at least. 
Now there exist a lot of conclusions on Wallis ratio, its variants, generalizations 
and applications, for example, [2, 11, 21, 23, 24, 28, 30, 31, 32, 33, 34, 35, 40, 41, 
44, 47, 48] and related references therein. 

In [26], D. Kershaw proved a double inequality 



s 

^+2 



r(x + i) 



r(x + s) ^ 

for < s < 1 and x > 1. It is clear that the inequality (5) can be rearran 

"r(x + i)T'/('-^) 



(5) 



as 



< 



V{x s) 



-X < 



This suggests us to introduce a function 



r{x + 1) 



-I l/(t-s) 



T{x + s) 



s 1 
s = t 



(6) 



(7) 



on X G (—a, oo) for real numbers s and t and a = min{s, t}. 

In [13, 22, 34, 35, 36, 39, 51, 52], the monotonic and convex properties of Zs,t{x) 
were established by using Laplace transform and other complicated techniques. 
Their basic calculation is as follows: 



z's,t(.x) = [Zs,t{x)+X 



V;(,T + t) - ip{x + s) 



t-s 

2 

+ 



ip'{x + t)- iIj'{x + s) 



[i^ix + t)- ij{x + s)f + {t- s)[^p'{x + t)~ ^P'{x + s)] 



This further suggests us to consider the following two functions: 



ip^x + t) — -tpix + s) 
A,,t(a;) = { [ t^s 

[V''(x + s)]2+V'"(a: + s), 



ip'ix + t) --ip'ix + s) 
t - s ' 



s 1 
s = t 



and 



Qs,t{x) = bP{x + t)- tp{x + s)]2 + {t- s)[V''(a; + t)- + s)] 

on a; e {—a, oo) for real numbers s and t and a — min{s, t}. 

In [3, p. 208], [8, Lemma 1.1] and [9, Lemma 1.1], the inequality 

Ao,oW = [V'W]2 + V"W>0 



(8) 
(9) 

(10) 

(11) 
(12) 
(13) 
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on (0,00) was verified. In [4, 8, 9], this inequality was applied to provide some 
symmetrical upper and lower bounds for T{x) in terms of tp{x) as follows: 

exp{a[e'^(-)(V'(x) - 1) + l]} < ^ < exp{p[e^^^\i;{x) - 1) + l]}, (14) 

where x* = 1.4616 • • • denotes the only positive zero of tp{x), and a and /? are real 
constants. 

The first aim of this paper is to present the completely monotonic property of 
the functions As t{x) and 8,5 4(2;) on (—a, 00), which implies the positivity of the 

function Ao.o(a;) defined by (13). 

Theorem 1. The functions As,t{x) for \t — s\ < I and —As,t{x) for \t — s\ > 1 are 
completely monotonic on x G {—a, 00). So are the functions Qs,t{x) for \t — s\ <1 
and —@s,t{x) for \t—s\>l on x G {—a, 00). 

The second aim of this paper is, by making use of Theorem 1, to provide an 
alternative proof for the monotonicity and convexity of the function Zs,tix), which 
is quoted as follows. 

Theorem 2 ([13, 22, 36, 39, 51, 52]). The function Zs,t{x) in {—a, 00) is either 
convex and decreasing for \t — s\ < 1 or concave and increasing for \t — s\ > 1. 

It is well-known [55] that Wallis cosine or sine formula is 

f.7r/2 /■it/2 



r-TT/Z r-TT/A 

/ sin" xdx = / cos" x Ax 
Jo Jo 

(tt (n-1)!! 

_ 0Fr((n + l)/2) ^ 



nT{n/2) 



2 n!! 
(n- 1)!! 



n 



II 



for n even, 

(15) 

for n odd, 



where nil denotes the double factorial. It has been estimated by many mathemati- 
cians and a lot of inequalities were established in, for example, [12, 14, 15, 16, 17, 
18, 19, 20, 25, 27, 43, 54, 58] and related references therein. 

The third aim of this paper is. by utilizing Theorem 2, to prove two sharp 
double inequalities relating to Wallis cosine or sine formula (15) and to bound the 
probability integral or error function as follows. 

Theorem 3. For n S N, 

1 < (2^-1)" < 1 (16) 

V'7r(n + 4/7r - 1) " (2n)!! ^7r(n + l/4)' 

^ <J^t.„< (17) 



2^n + 97r/16- 1 (2n+l)!! 2^/n + 3/4 



and 



< / e-^ dx < ^ ^ (18) 



v/l + (97r/16-l)/n J-^ ^1 - 3/(4n) 

In particular, taking n ^ 00 in (18) leads to 



f 

J — C 



e"^ dx = s/^. (19) 



The constants ^ — 1 and \ in (16) and the constants ^ — 1 and | in (17) are the 
best possible. 
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Let 



nn = —, TT (20) 

r(l + n/2) ^ ' 

be the volume of the unit baU on M". The fourth aim of this paper is, by empfoying 
Theorem 2, to recover a double inequality for ratio of the volumes of the unit balls 
in R""-^ and M" respectively as follows. 

Theorem 4 ([2, Theorem 2]). For n e N, the inequality 



27r f7„ - V 27r ^ ' 

holds if and only if A < i and B > ■§ — 1 • 

The final aim of this paper is, by using Theorem 1, to generalize the inequal- 
ity (14) to a monotonic property as follows. 

Theorem 5. For real numbers s and t, a = min{s,t} and c S (— a,oo), let 

'T{u + t) r(c+s)" 



-f 



In 



_r(M + s) T{c + t) 

/ [ip{u + s) - V'(c + s)] du, 

^ J c 

on X G (— a,(X)). Then the function 



du, s ^ t 
s = t 

X ^ c 



(22) 



r , , J Kt(a;)-l]exp[5;t(a;)] + r 

fsA^) = \ ^ (23) 

on X G {—a, oo) is decreasing for |s — i| < 1 and increasing for |s — t| > 1. 

Remark 1. If taking c = x* , then the case s = t in Theorem 5 becomes [4, Theo- 
rem 4.3]: For < a < 6 < oo and x G (a, 6), the inequality (14) holds with the best 
possible constant factors 

^^{Q{h), b<oo ^^^^^^^ ^24) 

1, b= oo 



where 



lnr(x)-lnr(x*) 



Q(,)^^e'A(^)[V.(.)-l] + l' (2^^ 
2. Lemmas 

In order to prove our theorems, the following lemmas are necessary. 
Lemma 1. Let f{x) be defined in an infinite interval I . If 

lim f{x) = 5 and f{x) — f{x -f e) > 

for some given £ > 0, then f{x) > 5 on I. 
Proof. By induction, for any x £ I, we have 

f{x) > fix + e)> f{x + 2e)> ■■■> f{x + ke)^5 
as fc — > oo. The proof of Lemma 1 is complete. □ 
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Lemma 2 ([1]). For any positive integer n G N and x > 0, 

"fl 1 " 

'0 



r 

ij}{x) = Inx + / 
Jo 



u 1 — e 



(26) 



As X ^ oo, 



^(")(x) = (-ir+W -__e-™d«, (27) 
JO ^ 

V("-^)(a: + 1) = ^(^-^\x) + ~ (28) 

hix -- <ih{x) <\5ix (29) 

X ' 2x ^ ' 



Lemma 3 ([37, 38]). For s > r > 0, 

exp [{s - r)^(s)] > > exp [(s - r)V'(r)] . (31) 

Lemma 4 ([57]). ^ product of finite completely monotonic functions is also com- 
pletely monotonic. 

3. Proofs of theorems 
Now we are in a position to prove our theorems. 
Proof of Theorem 1. Direct computation and utilization of (28) gives 

Qs,t{x) - es,t{x + 1) = {[i}{x + t)+ll){x + t+ 1)] - [i}{x + S)+Il){x + S + 1)]} 

X {[i}{x + t)-i}{x + t + 1)] - [i){x + s)-i}{x + s + 1)]} 

+ (i - s){W{x + i) - V'(a; + t+l)]- [/(a; + s) - i}'{x + s + 1)]} 
f [^(a; + i + 1) + V(a; + f)] - [V^a; + s + 1) + V(a; + s)] 



2x + s + t \ {t-sf 



{x + s){x + t) J {x + s){x + t) 

{t-sf 



(32) 



As,t(a;) 



(a; + s)(a; + i) 



and 



hs,t{x) - As,t{x + 1) = — — + 



t—s\x+s x+s+1 x+t x+t+1 

+ 2{s + t + l)x + s'^ +t'^ +s + t 
{x + s){x + s + l){x + t){x + i + 1) 

^ 

{x + s)(a; + s + + i)(a; + t + 1) ' 

Since lima^^oo A^s\ {x) = for any nonnegative integer i by (26) and (27) and the 
function 

As,t(a:)-As,t(x + 1) 

1 - (S - f)2 

is completely monotonic by Lemma 4, that is, 

, [A,,t(a;)-A,,t(a; + l)]W _ (-l)'Ai^^(a;) - (-l)^Ai;^(a^ + 1) 
^ ^ l-(s-i)2 l-(s-i)2 
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(-l)'A*'' (x) 

on {—a,oo), then i_(s^'t)2 > follows from Lemma 1. This means the function 

i-''(s^ty^ is completely monotonic in (— q:,oo). 



Since the function 



(x+s)(x+t) 



is completely monotonic, then the function 



e,,t(a;)-e,,t(.T + l) 



is completely monotonic on (— a, oo) by considering (32) and Lemma 4, which is 
equivalent to 

"e,,t(x)-e,,t(x + i) 



(-1)* 



w (-i)feei^>(x)-(-i)'-ei^,^(:r + i) 

1 - (S - i)2 



> 



for nonnegative integer k. Further, from lima;^oo ©s t (x) = for nonnegative inte- 
ger k, which can be deduced by utilizing (26) and (27), and Lemma 1, it is concluded 



that 



i-l)''@i''}(x) 



> for any nonnegative integer k. This implies {—l)''@l'^}{x) ^ 



l-(s-t)2 

if and only if |i — s| ^ 1. Therefore, the functions Qg.tix) for — s| < 1 and 
—Qs,t{x) for |t — s| > 1 arc completely monotonic on (— a, oo). 

Since Qs,t{x) = {t — s)^As,t{x), the function As,t{x) has the same monotonicity 
property as Qs,t{x) on (— a,oo). The proof of Theorem 1 is complete. □ 

Proof of Theorem 2. It is clear from (9) and (10) that 



(33) 



for t ^ s. By Theorem 1, it is easy to see that 6s,t(a;) ^ and As^t(x) ^ in 
(— a, oo) if and only if |f — s| ^ 1. Then z'J^f{x) ^ for |f — s| ^ 1 follows from 
formula (33). The convexity and concavity of the function Zs,t{x) is proved. 
The inequality (31) is equivalent to 



max-^ e 



'}> 



r(r) 



> mm i e 



} 



for any positive numbers s > and t > 0. This implies 



tix) 



T{x + 1) 
T{x + s) 



l/{l-s) 



H-r + t) - uix + s) 



< e 



t-s 

^(x+t) ^{x + t)- jjjx + s) _ ^ 
t - s 

^(^+*V'(a;+<e)-i 



(34) 



< V'(a; + i)e^(^+*' - 1 



and 



z'sAx) > e 



+s)'^i^ + t)-'^ix + s) 



- 1 



= e'^(^+''V'(a; + - 1 
> V'(a; + s)e'^(^+") - 1 

if assuming t > s > without loss of generality, where ^ G (s, t). 
By the inequality (29), we obtain 

x^l)'{x)e-'^/'' < V''(a:)e'^(^) < a;^'(x)e-i/2^ 

for a; > 0. Using the asymptotic representation (30) yields 

lim \x%/j'{x)e-^/''] = 1 and lim \xi/j' {x)e-'^/^''] = 1. 



(35) 



(36) 
(37) 
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Hence, 



lim W{x)e'^^''^ = 1. (38) 
Combining (38) with (34) and (35) leads to 

lim 4 t{x) < lim W{x + i)e^(^+*)l - 1 = lim \ip'{x + t)e'''(^+*)] -1 = 



x+t—^oo 



and 



lim z', t{x) > lim fV/(x + s)e'''(^+")l - 1 = lim (x + s)e'^(^+^)] -1 = 0. 

X — *OC' ' X — >oo x+s — ^oo 

Thus, it is concluded that lima;^oo z'^ ^{x) = 0. 

Since z'^j.{x) ^ on x e (— a,oo) for \t — s\ ^ 1, then the function z'gj.{x) 
is increasing/decreasing on x S (— a,oo) for |f — s| ^ 1. Thus, it follows that 
z's^ti^) ^ and Zs,t{x) is decreasing/increasing in a; e (— a,oo) for |f — s| ^ 1. The 
monotonicity of the function Zsj,{x) is proved. □ 

The second proof of convexity of Zgj {x). It is sufficient to show the function 



i!{x - 



ipix + t) 



s-t 



T{x + s) 



T{x + 1) 



l/(s-t) 



s 1 



s = t 



(39) 



on {—a, oo) is increasing for |t — s| < 1 and decreasing for \t — s\ > 1. 
Straightforward calculation yields 

^^ip{x + s) - t{j{x + t) _^ In r(a; + s) - In r(a; + t) 



ln$s,t(x) = <^ s-t 

hiip'{x + s) + ip{x + s), 



t 



and 



[ln$,,t(x)]'= < 



+ 



1p' {X + S) - Tp' {X + t) ^ 1p{x + S) - ll^ix + t) 

tp(x + s) — ip{x + t) 
il)"{x + s 



= < 



il)'{x + s) 

s - t 



iiix + s) + i^ix + t) 



+ %l}'{x + s), 

As,t(a;), s^t 



s 1 
s = t 

ST^t 

s = t 



. %1!'{X + S) 

In virtue of Theorem 1 , it is concluded that 



t. 



[ln$.,t(a;)]' 



>0, if|i-s|<l, 
<0, ifk-s|>l. 



□ 



(40) 



The second proof of convexity of Zs_tix) is complete. 

Remark 2. In [6, 7, 8, 9, 22], the inequahties 

^'(a;)e'^(^) < 1, x>0 

and (13) were proved and used to construct many inequalities for bounding the 
gamma function T{x), the psi function ip and the trigamma function il>'{x) such 
as (14). 

In [22], as a corollary of Theorem 2, the following inequality was deduced: 



r(x + t) ]'/(*-^^ 

r(a; + s) 



< 



t-s 



V'(a; + t) — '>p{x + s) 



(41) 



holds for < |f — s| < 1 and with reversed sign if |t — s| > 1. 
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The second proof of convexity of -Zs,t(a;) implies that inequahties (13), (40) 
and (41) are equivalent to each other. 

Remark 3. It is conjectured that the function ^s,t{x) for |f — s| > 1 and its recipro- 
cal for |t — s| < 1 arc logarithmically completely monotonic on a; S {—ct, oo), which 
modified an open problem posed in [51, 52]. 



Proof of Theorem 3. Let the sequence 9i (n) be defined by 

(2n- 1)!! _ 1 
(2n)!! 

for n e N. In order to obtain the inequality (16), it is sufficient to show 

^ <e^{n)<--l 

4 TT 

for all n e N. Indeed, formula (15) and Theorem 2 implies 



ei{x) 



T{x + l) 



X = Zi/2,l{x) 



(42) 
(43) 
(44) 



_r(a; + l/2)_ 

is convex and decreasing in (— ^,oo), and then the sharp double inequality (43) 
can be deduced by observing that ^^i(l) = ^ — 1 and limr,.^^ Oi{x) = j. 
Let the sequence 62{n) be defined by 

(2n)!! 



(2n+l)!! 2Vn + ^2(n) 
for n e N. In order to obtain the inequality (17), it is sufficient to show 

for all n e N. Indeed, formula (15) and Theorem 2 implies 



02ix) = 



r(x + 3/2) 



1 2 



-x = ^1,3/2(2;) 



(45) 
(46) 
(47) 



r(a; + 1) 

is convex and decreasing in (—1, 00), and then the sharp double inequality (46) can 



be deduced by observing that 02 (1) 



16 



— 1 and lima; 



The rest is the same as the proof of [12, Theorem 1.2]. 
Proof of Theorem 4- The inequality (21) can be rearranged as 



.02{X) = 1 



A /n\ 
2 < 1 2 J 



< 



B 



□ 



(48) 



for n E N. Since 21,1/2(2;) is decreasing in (—5,1) by Theorem 2, considering 
^1,1/2 (^) = f ~ 1 liniai^oo -^1, 1/2(2;) = J, then the inequality (48) is concluded. 

□ 



Proof of Theorem 5. Straightforward computation yields 



1 



In 



T{x + t) r(c+s) 



r(x + .s) r(c + 1) 



^ip{x + s) — ipic + s), 
It is clear that gs^t (c) = (?s ^ (c) = and 

r 2l;{x + t)- 4){x + s) 

9'lt{x)=\ t-s 
[i)'{x + s) > 0, 

This leads to (a; — c)g'g ((a;) > and gs,t{x) > fov x ^ c. 



s = t. 



(49) 



> 0, Sjtt; 
s = t. 
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Differentiation yields 

{[<,(a.)-l]exp[ff;,(x)] + ip' 

where 

hs,tix) = WsA^) - 9s,t{x)g'^^t{x) - 1] e■y^^[g'^ ^{x)] + 1. 
Since hs,t{c) = 0, [g'^^t{x)f + g"'tix) > by Theorem 1 and 

KA^) = -9sAx){[g':A^) + g':'A^)]^}^M9's,ti^)] ^ o 

for |s — t| ^1 and x c, we obtain {x — c)hsAx) ^ and hence f'sAx) ^ for 
\s — t\ ^1 and X ^ c. □ 
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